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The ground states of Bose-Einstein condensates of spin-2 bosons are classified into three distinct 
(ferromagnetic, "antiferromagnetic", and cyclic) phases depending on the s-wave scattering lengths 
of binary collisions for total-spin 0, 2, and 4 channels. Many-body spin correlations and magnetic 
response of the condensate in each of these phases are studied in a mesoscopic regime, while low- 
lying excitation spectra are investigated in the thermodynamic regime. In the mesoscopic regime, 
where the system is so tightly confined that the spatial degrees of freedom are frozen, the exact, 
many-body ground state for each phase is found to be expressed in terms of the creation operators 
of pair or trio bosons having spin correlations. These pairwise and trio- wise units are shown to bring 
about some unique features of spin-2 BECs such as a huge jump in magnetization from minimum to 
maximum possible values and the robustness of the minimum-magnetization state against an applied 
magnetic field. In the thermodynamic regime, where the system is spatially uniform, low-lying 
excitation spectra in the presence of magnetic field are obtained analytically using the Bogoliubov 
approximation. In the ferromagnetic phase, the excitation spectrum consists of one Goldstone mode 
and four single-particle modes. In the antiferromagnetic phase, where spin-singlet "pairs" undergo 
Bose-Einstein condensation, the spectrum consists of two Goldstone modes and three massive ones, 
all of which become massless when magnetic field vanishes. In the cyclic phase, where boson "trios" 
condense into a spin-singlet state, the spectrum is characterized by two Goldstone modes, one single- 
particle mode having a magnetic-field-independent energy gap, and a gapless single-particle mode 
that becomes massless in the absence of magnetic field. 

PACS numbers: 12.20.Ds, 42.50.Ct, 42.50.Lc 



I. INTRODUCTION 
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Bose-Einstein condensates (BECs) of trapped alkali atoms have internal degrees of freedom due to the hyperfine 
of the atoms. When a BEC is trapped in a magnetic potential, these degrees of freedom are frozen and the state 
the BEC is described at a mean-field level by a scalar order parameter. When a BEC is trapped in an optical 
^^otential, however, the spin degrees of freedom are liberated, giving rise to a rich variety of phenomena such as spin 
■^domains |^ and textures Here the order parameter has 2f -\- 1 components that transform under spatial rotation 
C^s the spherical tensor of rank /, where hf is the hyperfine spin of bosons. 

^ ' Mean field theories (MFTs) of spinor BECs were put forth for both spin-1 |^-^ and spin-2 cases. According 
I to them, the m = magnetic sublevel of an antiferromagnetic BEC is not populated in the presence of magnetic field 

■ Ofor both spin-1 and spin-2 cases. However, Law et al. j|] used many-body theory to show that the m — sublevel of a 
^pin-1 BEC is macroscopically populated due to the formation of spin-singlet "pairs" of bosons. It was subsequently 
^hown that the m = sublevel of a spin-2 BEC is also macroscopically occupied due to the formation of spin- 
• «inglet "trios" of bosons. The physics common to both cases is that the spin-singlet state is isotropic and therefore 

. ^each magnetic sublevel shares the equal population. 

k> , A question then arises as to what extent and under what conditions MFTs are applicable. It is now understood |^,^ 
5_jthat the validity of MFTs is quickly restored with increasing an applied magnetic field. Thus for the many-body spin 
Cdcorrelations to manifest themselves, the external magnetic field have to be very low. 

The spin-singlet pairs of bosons should be distinguished from Cooper pairs of electrons or those of "^He due to the 
statistical difference of constituent particles. The Cooper pairs consist of fermions, so that the state is symmetric 
only under the permutations that do not break any pair. For the case of spin-singlet pairs of bosons, the state is 
symmetric under any permutation of the constituents. The Bose-Einstein statistics leads to a constructive interference 
among permuted terms, giving rise to a highly nonlinear magnetic response to be discussed later. In contrast with 
usual antiferromagnets, where antiparallel spins are alternately aligned (Neel order), "antiferromagnetic" BECs do 
not possess such a long-range spatial order because the system lacks crystal order. The antiferromagnetic phase of 
BECs is also called polar 0]. 
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In Refs. HJ^,^, only spin degrees of freedom are considered by assuming that the spatial degrees of freedom are 
frozen. In this paper, we relax this restriction and develop a theory of spin-2 BECs that enables us to study many- 
body ground states and the excitation spectrum thereof on an equal footing. For spin-1 BECs, this program has 
been carried out in Ref. |l(]f| . Many-body spin correlations and magnetic response of BECs, the results of which were 
briefly reported in Rcf. M, are expounded. The role of symmetry of the ground state in determining the character of 
the excitation spectrum is also elucidated. 

This paper is organized as follows. Section || derives an effective Hamiltonian that enables us to study many-body 
spin correlations and low-lying excitation spectrum of spin-2 BECs on an equal footing. Section [II reviews mean-field 
properties of each phase of BECs. Section studies many-body spin correlations and magnetic response of BEC. 
The energy eigenstate is explicitly constructed using the creation operators of boson pairs and trios. The degeneracy 
of the eigenstate is examined and some novel magnetic response such as a huge jump in magnetization and the 
robustness of the minimum-magnetization state against an applied magnetic field are discussed. Section ^ derives 
excitation spectra of BECs using the Bogoliubov approximation. All excitation spectra are obta ined analytically and 
the relations of their characters to the symmetry of the ground state are discussed. Section VI summarizes the main 
results of the present paper. Appendix |A| recapitulates the parametrization of the order parameter of spin-2 BECs, 
and appendix B| describes a method of calculating Zeeman-level populations. 



II. FORMULATION OF THE PROBLEM 



A. Interaction Hamiltonian 

Consider a system of identical bosons with hyperfine spin / and let ^^{r) (m = /, / — 1, • • • , — /) be the field 
operator that annihilates at position r a boson with magnetic quantum number m. The field operators are assumed 
to obey the canonical commutation relations 

[^„(r), (^')] - S^nrAr - r'), [*,„(r), #„(r')] = 0, Hir')] = 0, (1) 

where the Kronecker's delta Smn takes on the value of 1 if m = n and otherwise. The Bose-Einstein statistics 
requires that the total spin of any two bosons whose relative orbital angular momentum is zero be restricted to 
F = 2/, 2/ — 2, • • • , 0. We may therefore use F as an index for classifying binary interactions between identical 
bosons: 

J2 (2) 

F=0,2,---,2f 

where F'^^ describes an interaction between two bosons whose total spin is F. To construct V''^\ consider the 
operator AFMir, r') that annihilates at positions r and r' two bosons with total spin F and total magnetic quantum 
number M: 

f 

AFM{r,r')^ (F,M|/,mi;/,m2)*™,(r)*„,(r'), (3) 

mi ,m2 — — / 

where {F, M\f,m,i; f,m2) is the Clebsch-Gordan coefficient. We may use ApM to construct V'-^^ as 

p 

V^^^ = \(dr (dr'v^^\ry) ^ i^lr, r')iFM(r, r'), (4) 

where 7;'^^(r,r') describes the dependence of the interaction on the positions of the particles. Because of the com- 
pleteness relation X^fa/ l^i (^i -^'^I = Ij where 1 is the identity operator, we find that 

F 

E E Ay,,{ry)AFM{ry) = :n{r)h{r'):, (5) 

F=0,2,---,2f M=-F 

where 
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(6) 



is the total density operator and :: denotes normal ordering, that is, annihilation operators are placed to the right of 
creation operators. Integrating Eq. (^) over r, r' yields 



fdr f dr' ^ A^^^,{r,r')ApM{r,r') - 7V(7V- 1), 

•' F.M 



(7) 



where N = J n{r)dr is the total number of bosons. 

In the case of a dilute Bose-Einstein condensate of neutral atoms, we may to a good approximation assume that 
v^^\r, r') = gFS{r — r'), where gp characterizes the strength of the interaction between two bosons whose total spin 
is -F, and is related to the corresponding s-wave scattering length ap as 



9F 



-ap- 



Equation M) then becomes 



^ fi r IT" 



(8) 



(9) 



M=-F 



In the following discussions we shall focus on this case and therefore denote ApMir^r) simply as ApMir). 
When f = 2, F can take on values 0, 2, and 4. For F = 0, we have 



y{0) ^ 90 

2 



drAlAr)Aooir), 



where 



For F = 2, we have 



^o(^) 



V5 



2*2(r-)*-2(r) - 2*i(r)*_i(r) + ^{r] 



= | jdr A\^,{r)A,M{r). 



(10) 



(11) 



(12) 



M=-2 



For F = 4 we have 



V 



(4) _ 34 

2 



dr Al;,j(r)AiM{r) 

nT__4 

2 

n^{r):-Al,{r)Aoo{r)~ ^ Al,,{r)A2Mir) 



'4 Idr 



(13) 



where Eq. (||) was used in obtaining the second equality. Summing Eqs. (|lO|), (^ and (|T^), we obtain the interaction 
Hamiltonian as 



' dr 



34 : n^{r) : +(go - 54)^^0 (^)^oo(r) + [92 - 54) XI ^2M('^)^2Af (r) 

To eliminate the last term in Eq. (p[), we note the following operator identity: 

2 

i : F (r) : + ^ 4,,(r)i2M(r) + A^r) A,,{r) = ^ : ^^(r) :, 



(14) 



(15) 



M=-2 
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where F — {F^, Fy,F^) represents the spin density operators defined by 

1 

m,n— — 1 

with (i — X, y, z) being the (m, n) -components of spin-2 matrices given by 



(16) 



XX _ 

■' ~ 2 
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-2 J 



(17) 



We may use Eq. ( |T5| ) to eUminate the last term in Eq. (|TJ), obtaining 

V=^ jdr\co: n^{r) : +ci : F^(r) : +C2iSo(^)^o(r-) 
where co = (4^2 + 3g4)/7, ci = (^4 - 52)/7, and C2 = (7^0 - 10^2 + 354)/7. 



(18) 



B. Total Hamiltonian 

In the foUowing discussions we shall assume that the external magnetic field is weak enough to ignore the quadratic 
Zeeman effect. Then the total Hamiltonian H of the system consists of the kinetic energy term Hke, the trapping 
potential energy term i/pE, the linear Zeeman term TJlz, and the interaction term V: 



m = — 9 \ ^ 



m=-2 
2 



m=-2 

^LZ =-pJdr J2 *m/m„*n = -p J dr J2 "^^n 



(19) 
(20) 
(21) 



where p (> 0) is the product of the gyromagnetic ratio and the external magnetic field which is assumed to be applied 
in the z-direction. 

When the system is spatially uniform, i.e. Utia.p{f) = 0, it is convenient to expand the field operator in terms of 
plane waves: 



Then the single-particle part of the Hamiltonian becomes 

Hq = Hke + Hi^z = Y^'^''~ 



k.m 



where Ck = h^k^ /2M ^ and the interaction Hamiltonian becomes 



(22) 



(23) 



(24) 



where 
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Pk= I n{r)e *'° ''dr = ^ „,aq+fc,,„, (25) 

Fk = J F(r)e-*'= ''dr = ^ f™„a^,„aq+fc,n, (26) 

qmn 

ifc = y" ioo(r-)e'*'-dr = i= ^(-l)'"aq,™afc_q,_„. (27) 



III. MEAN-FIELD THEORY 



Mean field theory (MFT) of spin-2 BEC was discussed in Refs. [|6|,[7|. We present here a brief summary of as much 
of this theory as is relevant to later discussions. We shall also present some new results. When the system is uniform, 
BEC occurs in the fe = state. We therefore assume the following trial wave function 

N 



(28) 



where the complex amplitudes are assumed to satisfy the normalization condition 

2 

E ICml' = l- (29) 

m=-2 

The variational parameters Cm are determined so as to minimize the expectation value of H over the state (|2 

(H) = ^N{N - 1) + ^N{N - l){fy + ^N{N - l)\{s.)\' - pN{^), (30) 

where 

(/) — fmnCmCn, (31) 

mn 

(A)^EHCml', (32) 

m 

(s_)^iE(-l)"C™C-™. (33) 

m 

When the external magnetic field is applied in the z-direction, only the z-component of (/) is nonzero. We thus 
obtain 

{H)^^NiN^l) + ^^eo, (34) 

where 

eo = ci{h)' + %\{s-)\^-p{h) (35) 
5 

withp= 2Vp/{N- 1). 

The mean-field solution should be determined so as to minimize eo subject to the normalization condition (p9|): 

^ -\ (2ci(A)-p)mCm-AC™ + ^C2(-l)"(s_)Clm=0, (36) 



where A is a Lagrange multiplier. Multiplying both sides by (— l)™C_„i and summing over m yield 

A - ^02^ = 0. (37) 



5 

On the other hand, multiplying both sides of Eq. (^6|) by Cm and summing over m yields 



(2ci(/.) -p)(/.) - A + ^C2|(J_)P = 0. (38) 
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A. Ferromangetic BEC 



When (s_) = 0, Eq. ® gives 



(2ci(A) -p)m- A 



For nonzero components C„j 7^ 0, Eq. (39) gives 



(2ci(A) -p)m-A = 0. 



(39) 



(40) 



The case oi p — 2ci{fz) wiU be discussed in Sec. [II C. When p ^ 2ci(/z), only one component can be nonzero. The 
mean- field solutions, magnetizations, and mean- field energies are therefore given by 



C = 6^^(1,0,0,0,0), (A) =2, eF = 4ci-2p, 
C = 6^^(0,1,0,0,0), (/.) = 1, . 



6^^(0,0, 0,1,0), (/.) = -i, . 

C-e^*(0,0,0,0,l), (/.) = -2, eF = 4ci + 2p, 



(41) 
(42) 
(43) 
(44) 



where is an arbitrary global phase. 

The ground state is degenerate with respect to the global phase 
conservation of the total number of particle and leads to a massless Goldstone mode, as will be shown in Sec. p/B 



This represents the gauge invariance , i.e 



The conservation of the spin angular momentum does not lead to a new Goldstone mode because in ferromagnets 
all spins are aligned in the same direction and therefore the total spin angular momentum has the same piece of 
information as the total number of particles. 



B. Antiferromagnetic BEC 

The antiferromagnetic (or polar) phase of a BEC is defined as the one having nonzero spin-singlet pair amplitude, 



^s_; 7^ 0. When (s_) 7^ 0, Eq. ^ gives A = 2c2/5. Substituting this and Eq. (|38|) into Eq. (|3^) yields 

AF _ C2 V I 



^0 = 



:(/.)■ 



With A = 2C2/5, Eq. (|§) leads to 

(2ci(A)-p) 



2ciTO^ - -C2 I (/z) -pm 



(45) 



(46) 



When p 7^ 2ci{fz), the solutions of Eq. (|4^) is that only (C2, C-2) or (d, C_i) or Co is nonzero. Determining these 
values using conditions (fz) — X^m'^ICmP and Eq. (|29|), we obtain the mean-field solutions and the corresponding 
magnetizations as m 



c = 



l + %i,0,0,0,e'x./i 



^ = _i=e^0(^O,yi7c^,O,e'^iyr^C^,o), (/.) = 
C =6^^(0,0,1,0,0), (A)=0 



2p 


6^^- 


C2 




4C1-C2/5' 




5 


4ci - C2/5 


P 




C2 




2(ci -C2/5)' 


5 


4(ci - C2/5) 


0, 


_ C2 

5 ' 







(47) 

(48) 
(49) 

The mean-field solutions ( [17| ) and ( |48| ) are degenerate with respect to two continuous phase variables, that is, the 
global phase (j) and the relative phase Xm = 0-m — 0m (fn = 1,2) between the two nonvanishing amplitudes C±m- 
Corr espo nding to these two continuous degeneracies, we expect to have two Goldstone modes, as will be shown in 

Sec. |vq . 

When p = 2ci(/z), Eq. ( p8| ) with A = 2c2/5 gives |(s-)| = 1/2, which, together with (36), leads to Cm = 
e^"^" (— l)™C-m- Hence we have (fz) — 0. This is possible only when the external magnetic field is zero. The 
corresponding order parameter is given by 
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gi(<^2 — 0o) g«(0i — 0o) g— — 0o) g— i(02 — 0o) 

This solution is degenerate with respect to five continuous variables: one global gauge ((/)o), two relative gauges {(t>2 — <t>o 
and 01 — (po), and two variables 6 and 1}} that specify the amplitudes of the order parameter. As a consequence of 
these degeneracies, we expect to have five (three density-like and two spin-like) Goldstone modes, as will be shown in 
Sec. 

C. Cyclic BEC 

The remaining possibility is the case in which (s__) = and {fz) = p/2ci. This phase will be referred to as cyclic 
phase. The energy of this phase is given from Eq. (|3g) by 

e[r = -ci(A)^ (51) 

Let us now parameterize the order parameter of the cyclic phase as it will be needed to find the Bogoliubov 
spectrum. There are four equations (six real equations) that restrict the order parameter of this phase, that is, 

EK™I' = 1' (52) 

m 

\Y.^-l^C.^C-.^=0, (53) 



(/,)=^m|C™|2 = ^, (54) 

m 

(/+) = 2(aci + C1C-2) + x/6(crco + CoC-i) = o. (55) 



We use the representation ( AlOj ) of the order parameter derived in Appendix |^ to analyze the cyclic phase. This 



representation automatically satisfies the normalization condition (^2|). To meet the condition (p3|), we note that 

= ^TrM^ = - CiC-i + C2C-2 = ^ sinxcos(0 - 0). (56) 
The condition (|5^) therefore requires cither 

X = Qov(j)-6 = 7r/2. (57) 

On the other hand, Eq. ( p5|) becomes 

(/+) = 2sin(5cosx [cosV'sin(6' + 7r/6) — jsinV'sin(0 — 7r/6)] = 0, (58) 

whence we obtain 

(i) X = 7^/2, or (u) 5 = 0, or (in) ip = 7r/2 and 9 = 7r/6, or (iv) -0 = and 6* = -7r/6. (59) 
From conditions ( |57| ) and (|5^), we find the following three solutions and the corresponding magnetization: 
C±2 = i(cos 6* ± cosx — zsin0sinx)e"^, C±i = 0, = -^(sin0 -f i cos0sinx)e"^, (/z) = 2cosxcos0, (60) 

C±2 = 5(^±cos(5cosx-§sinx)e^'^, C±i = ife'"^ sin ,5cos x, Co = 75 (5 + sin x) e*^, (/,) = V3cos<5cosx, (61) 
C±2 = i(^±cos(5cosx+tsinx)e'*, C±i = ^e'^siufJcosx, Co - (-^ + sinx) e^^, (/,) = VScosJcosx, (62) 



where a global phase 0, which is chosen to be a particular value in Eq.(A6), is recovered. While these solutions include 
three parameters, the condition (|^ leaves only two parameters free. It should be noted that these three solutions 
give the same ground-state energy, and hence are equally possible unless magnetization exceeds \/3- When it exceeds 
■\/3, only solution (|60|) is possible. 
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A remark is here in order. As the representation (AlO) is obtained by assuming the full isotropy of space, it does 
not cover the whole order parameter space in the presence of magnetic field. We should therefore keep in mind that 
the solutions (|60|)-(|6^) do not necessarily exhaust the whole order parameter space of the cyclic phase. 

In the following discussions we shall focus on the solution (BOh . Making the absolute square of each amplitude yields 



IC: 



±21 



-(l±(^)/2)^ |CoP = :,(l-(/.)V4). 



(63) 



Hence we obtain 



,±2 




where 



— sin X sin 9 
tan(p±2 = T. , tan 00 = sin^cotl 



cos 9 ± cos X 
It follows from this or by direct calculation that 

(t>2 + <l>-2 — 200 = ±7r 



(64) 



(65) 



(66) 



Because of this restriction, the ground state of the cyclic phase is degenerate with respect to at least two continuous 
phase variables. We therefore expect to have at least two Goldstone modes, as will be shown in Sec.VD. 



D. The phase boundaries 

In the absence of external magnetic field, the ground-state energies for the three phases are given from Eqs. (|l]), 
(g|), and (13) by 

e^=4ci, er = |, e^ = 0. (67) 

It follows that each phase is specified by 



ferromagnetic ci < and ci — C2/2O < 0, (68) 
antiferromagnetic C2 < and ci — C2/2O > 0, (69) 
cychc Ci > and C2 > 0. (70) 

In the presence of external magnetic field we define each phase as follows: 

ferromagnetic {f^) ~ 2, 
antiferromagnetic (s_) 7^ 0, 

cyclic (.s_) = and {f^) < 2. 

By directly comparing the energies in Eqs. (^), (^), and (|5|), we find that each phase is specified by the following 
conditions: 

ferromagnetic ci < p/4 and ci — C2/2O < p/4, (71) 
antiferromagnetic C2 < and ci — C2/2O > p/4, (72) 
cyclic ci > p/4 and C2 > 0. (73) 
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IV. MANY-BODY SPIN CORRELATIONS AND MAGNETIC RESPONSE 



In this section we study the case in which the system is so tightly confined that the coordinate part of the ground- 
state wave function 4>o{r) is independent of the spin state and solely determined by /^ke, ^^pe, and the spin- 
independent part of V; that is, (poir) is the solution (/) to the equation 



2M 



trap 



with the lowest eigenvalue e = eq- This assumption is justified if the second lowest eigenvalue ei satisfies 

ei-eo» bUci|7V/V^'=«,|c2|iV/F'=ff, 



(74) 



(75) 



where = (/ dr\(j)fif)-^ is an effective volume which coincides with V in Eq. ( p2[ ) for the spatially uniform case 
(i.e., C/trap = 0). When the condition ( [TSf ) is met, the field operator may be approximated as — o.m(f'o, where 
am is the annihilation operator of the bosons that are specified by the spin component m and by the coordinate wave 
function (j>o- The spin-dependent part of the Hamiltonian can then be written as 



where 



H = 



2C2 



S+S- - pTz, 



a„,a- 



ma'a„ 



(76) 



(77) 



A. Spectrum and degeneracy 

We first make some remarks on the properties of the operators S- = S]_ = (ao)^/2 — aia_i -I- 0,20,^2- The operator 
iS+, when applied to the vacuum, creates a pair of bosons in the spin-singlet state. This pair, however, should not be 
regarded as a single composite boson because S+ does not satisfy the commutation relations of bosons. The operator 
5+ instead satisfies the SU{1, 1) commutation relations together with Sz = (27V -(- 5)/4, namely, 

[Sz,S±\ = ±S±, = -2Sz, (78) 

where the minus sign in the last equation is the only distinction from the usual spin commutation relations. As a 
consequence, the Casimir operator 5^ that commutes with S± and Sz is given by 

S"^ = -S+S-+Sl-Sz. (79) 

Consider an eigenspace Ti^ of 5^ with an eigenvalue v. The requirement that Sj^S- ^ S^_—Sz~S^ must be positive 
semidefinite means that, in Jiu, the eigenvalues of the operator Sz has a minimum value (2iVo + 5)/4, where A^o is a 
nonnegative integer. For a state |0) that belongs to the minimum eigenvalue, the norm of must vanish; hence 

V = S{S — 1) with S = (2iVo -|- 5)/4. We thus obtain the allowed combinations of eigenvalues {S{S — 1), S'2} for 5^ 
and Sz such that 

5= (27Vo + 5)/4 (TVo = 0,1,2,...) (80) 

and 

Sz = S + Ns (TVs =0,1,2,...). (81) 
Here we have introduced quantum numbers iVs and A'o, where the operator 5+ raises TVs by one and the relation 

iV = 2iVs + ^0 (82) 
holds. We may thus interpret No, as the number of spin-singlet 'pairs', and A^o as that of the remaining bosons. 
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Exact energy eigenvalues of Hamiltonian ([76|) can be obtained as follows. The operators S± are invariant under 
any rotation of the system, namely, they commute with the total spin operator JF. The energy eigenstates can thus 
be classified according to quantum numbers Nq and Ng, total spin F, and magnetic quantum number Fz- We thus 
denote the eigenstates as \Nq, Ns, F, F^; A), where A = 1, 2, ... , gNo,F labels orthonormal degenerate states, that is, 

(A^o, Ns, F, Fz- A'lTVo, Ns, F, F^^X) = 5x',x. (83) 

The number of degenerate states gNo,F for a given set of {A'o, Ns, F, F^} will be referred to as the size of the eigenspaces 
for {Nq, Ns, F, Fz}. It will be shown to be independent of Ns and Fz below. The energy eigenvalue for the state 
\No,Ns,F,Fz;\) is given by 

^ = ^[^(^ + 1) - 6iV] + ^NsiN + No + 3)~ pFz, (84) 

where the relation 2Ns + Nq = N is used. 

The degeneracy gNo,F can be calculated as follows. First we show that gNa.F is independent of Ns and Fz. This is 
seen by the following relations, 

{No,Ns,F,Fz;>:\J'±j'^\No,Ns,F,Fz;\) = {F ± Fz){F ^ Fz + l)5x' ,x, (85) 

where J-± = J-x ± iTy , 

{No,Ns,F,Fz;X'\S+S^\No,Ns,F,Fz;\) = Ns{No + Ns + 3/2)Sy^^, (86) 

and 

(A^o, Ns, F, Fz; X'\S-S+\Nq, Ns, F, Fz; A) = {Ns + l){No + Ns + 5/2)Sy,x. (87) 

These relations implies that the sizes of the eigenspaces for {A'o, Ns ± 1, F, Fz ± 1} are not smaller than the size of 
the eigenspace for {A^'o, Ns, F,Fz}. The degeneracy thus depends only on A'o and F. Next, we introduce a generating 
function of gNo,F defined by 

oo oo 

G{x,y)^Y.T.9No,FX^''y''- (88) 

Aro=0 -F=0 

Let hpj p^ be the total number of states with a fixed number of bosons N and a fixed magnetic quantum number 
Fz. This is given by the total number of combinations of nonncgative integers {n_2, n_i, rig, ni, 712} that satisfy 
n_2 + n^i + riQ + ni + n2 = N and — 2n_2 — "--i + ^^l + 2712 = Fz. It follows that 

00 00 2 2 

JV=OF^=-oo {uj} j = -2 j = -2 

where we assume \y\ < 1 and \z\ < 1 to ensure the convergence of the series. Let hN,F,F^ be the total number of states 
for given N, F, and Fz. Because /itv.f.f^ is independent of the value of Fz, we shall denote it simply as h^^F- 
The quantity hN,F^ is written in terms of the sum of hj^^p as 

hN,F, = E ^N,F, (90) 
F{>\F,\) 

and hence /ijv,F = hpf p — his[ p+i. Let us extend the definition of /ijv.F to the negative values of F through this 
relation. It follows then from Eq. ( 89 ) that 

oo oo 2 

^ y: h^.Fz''y^^+^+' = (y - 1) n (1 - ^y'^'y"- oi) 

Af=OF=-oo J = -2 

The right hand side of this equation can be written as the sum of two fractions Gi{z, y) + G2{z, y), where 

^ , X yjl - zy^ + z^y^) 

"^^^^'^^ (l-Z2/3)(l_^y4)(i_^2y4)(i_^3y6) > 
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and 



G2{z,y) 



2^,2 



1 — zy + z y 
(l-z)(l-zZ/)(l-zV)(l-^V)' 



(93) 



Making Maclaurin expansions of Gi and G2 around z = y = and regrouping them in terms of the form 
we find that Gi consists only of the terms with F > 0, and G2 of those with F < 0. We thus obtain 



The quantity h^^p is written by the sum of the degeneracy gNo,F as 

hN,F = ^ gN-2Ns,F, 
Ns{<N/2) 



(94) 



(95) 



and hence we can write gNo.F = hj^r^^p — hNg^2,F, where we assume that h^i p ~ /i-2._F = 0. It follows then from 
Eq. dgl) that 



00 00 



Y: 9n.,fzW''^'^' = (1 - z^y')G,{z,yl 

No=0 F=-oo 

and we finally obtain an explicit form of the generating function G{x,y) defined by Eq. ( ^8| ) as 

1 — xy + x^y^ 



00 C30 

Y E 5^o.Fa:^"2/^ = y-^(l - x')G^{xy-\ y) = 

No=0 F=0 



{l-xy)il-xy^)il^x^)' 



(96) 



(97) 



The total spin F can, in general, take integer values in the range < F < 2No. However, from Eq. (|97|) we find 
that there are some forbidden values. That is, F — l,2,5,2A'o — 1 are not allowed when A^o — 3fc(fc G Z), and 
F — 0, 1, 3, 2jVo — 1 are forbidden when A^q = 3fc ± 1. An easier way to find the forbidden values is discussed at the 



end of Sec. IVB 



B. Energy eigenstates 

The energy eigenstates \No, Ng, F, F^; X) can be constructed using one-, two- and three-boson creation operators. 

Let us define the operator ^j"'^ such that it creates n bosons in the state with total spin F = f and magnetic 
quantum number F^ — f when applied to the vacuum. Such states are unique when n < 3. Among possible operators 

^y"'*^, we choose the following five operators for constructing the eigenstates: 

4'^^ = a2^ (98) 
4')^ = -^H? 2alali + 2alaU = ^1^+ (99) 

A "^ ' 



2^^ = ^[2V2a|4 - Vi{a\f] (100) 



1,^'^ = -^i=[\/2(4)^ - 3V2ai t + i^{a\fa^-2 + 3V34(aLi)^ - Qy/2a\ala}_^] (101) 

if ^ = -^[{a\f - V64aI4 + 2{a\faU. (102) 

Note that if ^ and if ^ do not exist because of the Bose symmetry. Note also that the operators i^"''^ commute 
with F|_. 

Consider a set B of unnormalized states, 

Ini2,n20,n22,n30,n33) ^ (a^)"- (if ^)"- (if ^)"- (if ^)"- (if ^)"- |vac) (103) 
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with ni2, n2o, n22, ^^3o = 0, 1, 2, . . . , oo and 7133 = 0, 1. The state |ni2, n2o, ^22, '^307 '^■33) has the total number of bosons 
N = ni2 + 2(7120 + ^22) + 3(n3o + 7133) and the total spin F = = 2(ni2 + 7122) + 87733. When N and F are given, 
7733 is uniquely determined through the parity of F, namely, 

7733 = F mod 2. (104) 

If we introduce the following two parameters 

^ = 27120 + 37730, (105) 

I/ = 7720+«30, (106) 

7720 and 7130 are uniquely specified by them, and the remaining 7712 and 7122 are also determined as 

7712 - F - iV + /i, (107) 

7722 =A^-^-M-^^33. (108) 

The parameter set {7112, 7720, ?^22, f^30j "■33} is thus uniquely specified by the set {N, F, ^^v}. Let us consider the two 
commuting observables 

77^ = 2al2«-2 + 2(^-1^-1 " "33) (109) 

= a^_^a^2 + ^{aLiCL-i - 7133) (110) 

with 77,33 = P mod 2. Let P{n^^n^) be the projection operator onto the simultaneous eigenspace of 77^ and 77^ 
corresponding to eigenvalues and 77^ , res pectively. Noting that Aq includes the term 020-2 and Aq includes 
a2(a-i)^, it can be seen from Eqs. (|9^)-(10S ) that 

-P(f^M,f^l/)|'^l2,'^20,f^22,?^30,'^33) 7^ if 77^ = ^ and 71^ = i/, (111) 

and also 

P{n^^, 7i;.)|77i2, 7720, ?^22, «30, "-33) =0 if > fi 01 > v. (112) 

If we order the pair {iJ-yiy) in the lexicographic way, namely, {^^v) < {^',1^') if fi < ^' or (/i = /i' and v < v'), the 
above equations imply that a state in B specified by {N, F, fi, u} is linearly independent of the set of states specified 
by {N, F, fi', v'} with (^', v') < {fi, v). The set B is thus a hnearly independent set of states. Let h'j^ p be the total 
number of states in B with N bosons, total spin F, and magnetic quantum number F^ = F. A generating function 
of h'j^ p is calculated as 

00 00 00 1 

N=0 F=0 ni2,"20,"22,n30=0 n33=0 

^^"^'^^ -y-'G^{xy-',y), (113) 



(1 - xy2)(l - x^){l - a;2?/2)(i _ 



where Gi is defined by Eq. (|9^). Compared to Eq. (|94|), we have h'j^ p — h^^p. This imphes that B is complete, 
namely, the set B forms a nonorthogonal basis of the subspace T-L^^p^^p) in which magnetic quantum number F^ is 
equal to total spin F. 

The energy eigenstates can be obtained by partially applying the method of Schmidt's orthogonalization to the 
nonorthogonal basis B. Here by "partially" we mean that eigenstates corresponding to different energies are or- 
thogonal, but that those corresponding to the same energy are not always so. Let us consider a series of subspaces 
Ti(^p,=p^ ~ Ti^"-* D Ti*-^-* D • • •, where Ti.'-^^^ is spanned by the states with quantum number Ng [see Eq. (81)] sat- 
isfying iVs > j. Let p(^' be the projection operator onto Ti^-^\ From the relation S+P'--^^ — p(^+^)5+, we have 
pa)(i(2)tp- ^ (i(j2)tpp(o)^ g-j^^g p(o)|0) = for any state |(/)) e B, we have 

(p(o) _p(„2o))|„^2^„2o, 7722, 7730, 7733) = 0, (114) 
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implying that Ti^^^ is spanned by all the states {|ni2, ^20, '^^22, ?^30; '^33)} satisfying 7120 > j- We can then construct a 
new basis B' made up of the states of the form 

where P(7Vg=o) = P'"^^ ~ -P''^'' is the projection onto the subspace with A^s = (the kernel of S-). It is easy to see 

* „ » -> 2 - 

that the states belonging to B' are simultaneous eigenstates of {N, S+S-,J^ , J^z}, and hence energy eigenstates. The 
energy eigenstates with Fz < F can be constructed by applying {^-)^~^' to the states of B' . 
To summarize, the energy eigenstates can be represented as 

(^_)^^(4^)^)"-P(^3.o)(4)"^M4''^)"^M4''^)"™(4'^^)"^1vac), (116) 

with ni2,n2o, 7122^^130 — 0, 1, 2, . . . , 00, 71,33 = 0,1, and AF — 0,1, . . . ,2F. These parameters are related to 
{No,Ns,F,Fz} as 

A^o = ?^l2 + 2n22 + 3^30 + 37133 (117) 

Ns = 7720 (118) 
= 27712 + 27722 + 37733 (119) 

Fz = F - AF, (120) 



and the corresponding eigenenergy is given by Eq. (g4|) . Note that the states defined in ( |116D are unnormalized, and 
the states having the same energy (i.e., those belonging to the same set of parameter values {Nq, Ns, F, Fz}) are 
nonorthogonal. 

The representation (|116|) of the energy eigenstates utilizes the operator Ai"-*^ that creates correlated 77 bosons having 



total spin /. It might be tempting to envisage a physical picture that the system is, like in ^He, made up of 77„/ 
composite bosons whose creation operator is given by ^j"^^- However, this picture is oversimplified. First of all, the 

operator Aj"'' does not obey the boson commutation relation. In addition, the projection operator P(jVs=o) in (116) 
imposes many-body spin correlations such that the spin correlation between any two bosons must have vanishing 
spin-singlet component. Note that two bosons with independently fluctuating spins have a nonzero overlap with the 
spin-singlet state in general. The many-body spin correlations of the energy eigenstates are thus far more complicated 
than what an intuitive picture of composite bosons suggests. On the other hand, as long as quantities such as the 
number of bosons, magnetization, and energy are concerned, the above simplified picture is quite helpful. By way of 
illustration, we provide an alternative explanation for the existence of forbidden values for the total spin F, which were 
found earlier using the generating function (p7|). For example, to construct a state with F = or = 3, composite 
particles with total spin 2 must be avoided, namely, 7712 = 7722 = 0. Then we have iVo — 3(7130 -|- 7733), implying that 
= or F = 3 is only possible when A^o = 3fc(fc e Z). For a state with F = 2 or F = 5, we have 7712 -I- 7122 = 1 
and A'o = 1-1- 7122 + 3(7130 -I- 7733) implying that iVo 7^ 3fc(fc G Z). The above simplified picture is also helpful when we 
consider the magnetic response as discussed below. 



C. Magnetic response 



We consider here how the ground state and the magnetization Fz respond to the applied magnetic field p. From 
Eq. (|^, we see that the minimum energy states always satisfy Fz — F when p > 0. The problem thus reduces to 
minimizing the function 

EiFz,N,) = ^[Fz{Fz + 1) - 6N] + ^Ns{2N - 2Ns + 3) - pF,. (121) 
For this purpose, it is convenient to consider the cases C2 > and C2 < separately. 



1. Case o/c2 > 



Let us first consider the ca se C2 > 0. 

When Cl < 0, the energy (121) is minimized when A^s = 0, A^o = N, and F ^ Fz = 2N, and the ground state is 



given by (a2)^|vac), that is, the system is ferromagnetic. This result is consistent with that obtained from MET. 
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When ci > 0, let us rewrite the energy as 



E{F,,Ns) 



Cl 



2V 



off 



pV 



cff 



Cl 



2C2 



TVs 



2iV + 3 



const. 



(122) 



The energy is thus lower when is closer to pV^^ /ci — 1/2 and when No, is smaller. In most of the parameter 
space, the ground state is |7Vg = 0, A^o = F, F^ = F; A) with F taking the allowed integer closest to pV°^ /ci — 1/2. 
Since these states satisfy {S+S-) = 0, they belong to the ferromagnetic phase or the cyclic phase. The separatrix 
between the two phases is given by pV°^ /ci — 1/2 ~ 2N — 1. We thus find that the ground state is ferromagnetic if 
Cl < / {2 N — 1 /2) ~ pV°^ / (2N ) and cyclic otherwise. This classification is consistent with the mean-field analysis 
given in Sec. [II D| . As seen in Sec. IV A , the above ground state is highly degenerated in general; this may originate 
from the continuous sy mmetr y that leaves free at least two paramet ers ch aracterizing the order parameter of the cyclic 
phase as shown in Sec. IIIC . According to the discussions in Sec. [V A , the degeneracy of the states \Ns = 0,Nq = 
F, Fz = F; A) is equal to the number of the combinations of {^12, n22, ^30} satisfying ni2 + 2n22 + 3n3o — N — 8^33 
and 2ni2 + 2^22 = F — 8^33. The number of trios, 7130, can take values in the range (TV — F)/3 < n^o ^ {N — F/2) /8. 
When the magnetic field is nearly zero and F ^ 0, there is little degeneracy and n^Q ~ N/3, namely, almost all bosons 
form trios. 

When pV°^ /ci — 1/2 is close to the forbidden values of (F = 2, 5 if A^o = (mod 8), and F = 0, 8 otherwise) for 
the above state with A^s — 0, F^ may take those values at the cost of increasing iVg to 1 or 2, since any of the three 
values (0,1,2) of iVo mod 3 is realized by setting Ng as 0, 1, or 2 [recall the relation = 2A's + Aq]- Whether or not 
the states |A^s = 1, A^o = A^ - 2, F, = F; A) and | A^s = 2, Aq = A^ - 4, F, F^ = F; A) can be the lowest-energy state 
depends on the ratio Ci/c\. In Fig |^, we give diagrams of the ground states for small pV"^ /ci. 



FIG. 1. Diagrams of the ground-state magnetization for C2 > and ci > 0. 



2. Case of C2 < 



In this case, it is convenient to introduce a new parameter 



20 

and consider the energy as a function of F^ and I = 2Nq — F^ 

( C2 \ 1 



Cl - (123) 



2V 



off 



F, 



off 



H 2 

-^^KZ + 2F + 6) + const. (124) 



Since C2 < 0, we see that E{Fz, V) is an increasing fimction of Z, namely, 

E{F,,l) > E{F,,l') for I > I'. (125) 

Let us consider the cases c'l < and c[ > separately. 

(a) c'l < — In this parameter region, MFT predicts that the system is ferromagnetic, namely, the lowest-energy 
state always shows (F^) = 2A^ regardless of the magnitude of magnetic field p. In the exact solution considered here, 
the magnetic response is different because of the offset term C2/ (SV^^) in Eq. ( |124| ). Since C2 < 0, this term counteracts 
the applied magnetic field. It is thus expected that magnetization is suppressed when magnetic field is weak. The 
exact ground state is derived as follows. When A^ is even, the state (i[,^^^)^/2|vac) has energy E{Fz = 0, Z = 0), and 
the state (a2)^|vac) has energy E{Fz = 2N,l = 0). Any other set {F^,!} gives an energy higher than one of these 
states. {A';^^^)^/^\va.c) is thus the ground state when E{Fz =0,1 = 0) < E{Fz = 2N, I = 0), or equivalently, 

pV'''' <c[(n+1]~^ (126) 



2^ 

and otherwise the ground state is (a2)^|vac) [see Fig |[. When A^ is odd, Fj = is attained only when Aq > 8, 
Fz = 1 is forbidden, and the state 02(^0^'^)^^""^^''^ jvac) has energy F(Fz = 2,1 = 0). It is easy to confirm that 
F(Fz = 2,1 = 0) < E{Fz = 0,1 = 6) always holds. Therefore, al{A''^^^)^'^-^'^^'^\vac) is the ground state when 
F(Fz = 2,1 = 0) < E{Fz =2N,l = 0), or equivalently. 



14 



pV'^^<c[^N+iy'f, (127) 

and otherwise the ground state is (a2)^|vac) [see Fig |[. These results indicate that in the parameter region of 
C2 ^ 8Nc'i < 0, magnetization of the ground state jumps from or 2 to 2N. Such a huge discontinuity does not 
appear in MFT with a hnear Zeeman effect. (However, in the presence of a quadratic Zeeman effect, such a jump 
occurs also in MFT [||.) 

FIG. 2. Diagram of the ground-state magnetization for C2 < and c'l < 0. 

(b) c'l > — Given Fz > 6, the minimum allowed value of I is determined as follows. Note that I = 2N — Fz — 
ANs is minimized when the number of singlets Ns is maximized. For Fz — 2N — 4fc (fc is an integer), the state 

(^Q^''^)'^(a2)^~^''|vac) gives 1 = 0. To increase Fz by one {Fz = 2N — 4fc + 1), one singlet pair must be broken and 

the minimum of Hs Z = 3 given by the state (A'^''^ )''~^ (a\)^^^''~^ Af"^^ \veLc) . Keeping the singlet part, Fz can be 

further increased to Fz = 2iV - 4/fc + 2 by the state (i[,^^^)'=-i(4)^"^''i2^^Vac) with I = 2. Since Fz = F = 2No-l 
is forbidden, Fz = 2N — Ak + 3 requires one more singlet pair to break up, resulting in I ^ 5 with the state 
^^(^2)t^fe-2(^t)A^-2fe-l^(2)t^(3)t|y^P^J^ When {pV^ + C2/8)/c[ - 1/2 falls between 2N - Ak and 2N - 4(fc - 1), the 
lowest energy is the minimum of e[Fz = 2N - Ak,l = 0), E{Fz = 2N - Ak + 11 = 3), E{Fz = 2N - Ak + 2,1 = 3), 



E{Fz = 2N - Ak + 3,1 = 5), and E{Fz = 2N - Ak + A,l = 0). From Eq. ( |124D , we expect that when \c2\/c[ is 
large, nonzero I pushes up the energy significantly and cannot be the ground state, so that Fz increases stepwise with 
AFz = 4. When \c2\/c[ is sm all, Fz will increase with the step size of AFz = 1. This is confirmed by explicitly 
calculating E{Fz, I) using Eq. (124), and we obtain the diagrams in Fig. ^. 



FIG. 3. Diagram of the ground-state magnetization for C2 < and c'l > 0. 

In the region pV^^ /c'l < fi{\c2\/c'i) with fi{x) = {x — 20) (9a; — 4)(32a;)^^ (see broken curves in Fig. ||), Fz increases 
by taking every integer. When pV^ /c[ > fi{\c2\/c'i), the values of Fz = 2N — 4A: -I- 3 are suppressed. In the region 
pV°^/c[ > /2(|c2|/ci) with f2ix) = (2:-20)(132:-20)(80x)-\ the values Fz =2N-Ak+1 are suppressed, and when 
pV^/c'i > /3(|c2|/ci) with f3{x) ^ {x- 20){x - 8){8x)-\ the values = 2iV - 4fc + 2 are further suppressed and 
the step size becomes 4. 



While the averaged slope AFz/Ap - V°^/c[ coincides with that in MFT, the offset term 02/(8^°*^) in Eq. (yj 
(see also the broken lines in Fig. ^ makes a qualitative distinction from MFT, namely, the onset of the magnetization 
displaces from p = to p = |c2|/(8y^). Note that the slope V^/c'i and the offset |c2|/(8y®'^) are determined by 
independent parameters. A typical behavior of the magnetic response when |c2| ^ c[ is shown in Fig. ^. 

FIG. 4. Typical dependence of the ground-state magnetization on the applied magnetic-field strength, for C2 < and c'l > 0. 



D. Property of ground states for C2 < 

We now calculate the Zeeman- level populations of the ground states for C2 < 0. In MFT, the lowest-energy states for 
C2 < have vanishing population in the m = 0, ±1 levels. In contrast, the exact ground states derived in the preceding 
subsection, (A[|^^^)^s (aj)"!^ (A^^^'^)"22 (^(''')t^»33 |vac) with 7122 = 0, 1 and 7133 = 0, 1, have nonzero populations in the 
m = 0, ±1 levels. The exact forms for the averaged population (aj,ja,„) are calculated as follows. The above ground 
states have the form of (^q^^^)"^^ |^) (x (iS+)^^|^) with |0) being a state with a fixed number (s = rti2 + 2n22 +37133) 
of bosons satisfying S-\(j)) = 0. The average Zeeman population for the ground states, 

^ ^ - (0|(5_)^sat^a,„(5+)^s|0) 

{aLo-m) = 7 (128) 

(0|(5_)^s(5+)^s|0) 

is then simply related to the average Zeeman populations for the state |0) as 

(aLa™) = (aLam)o H ^^^^{{(>La„^)o + (aL,„a-„i)o + 1), (129) 

s + 5/2 

where {al^am)o = {4>\0'm^m\4') / {4'\4>) ■ The derivation of the formula ( |129| ) is given in Appendix H. The formula 
imphes that when ATs 3> s, the Zeeman populations of the ground states is sensitive to the form of f^. 
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With this formula, it is a straightforward task to calculate average Zeeman-level populations for four types of 
ground states, (^t )ni2 (^(2)t)n22 (^^3)t)«33 i^a,,) with 7122 = 0, 1 and 7133 = 0, 1. The exact result will be given 

in Appendix A striking feature appears in the leading terms under the condition 1 <C ni2 <C A''s. The results are 
summarized as 

{a\ai) - (a^a_i) - iVs(l + ^33)7^12 (130) 

and 

(alao) ^ Ns{l + 2n22)/ni2. (131) 

As seen in the preceding subsection, with the increase of magnetic field, the ground state alternates among the four 
types of states. While this change causes a very small difference in magnetization, it leads to large changes in the 
average Zeeman-level populations, by a factor of 2 or 3. The origin of this drastic change may be explained as follows. 
Let us first consider the state (a2)"^^|vac) with ni2 > 0. This state has no population in Zeeman levels m = 0,±1. 
When A'^'^^ is applied to this state, the operator a\ that appears in ^2^^'^ has effectively a large amplitude of the 
order of ^7121. Hence the term 0.200 is dominant, and it approximately adds one boson to the m = 2 level and one 
boson to the m = level. Hence the m — population of the state A^2'^\a\)"'^^ |vac) is close to unity. This change is 
then amplified by a factor of Ns/ni2 according to the formula (129), leading to Eq. (131). Similarly, applying ^3 



effectively results in adding of two bosons to the m — 2 level and one boson to the m = —1 level through the dominant 
term {a\)'^a}_^. 

V. LOW-LYING EXCITATION SPECTRA 

In this section we study the low-lying excitation spectrum of spin-2 BECs in the thermodynamic limit using the 



Bogoliubov approximation. We shall see that the symmetry of each ground state discussed in Sec. Ill is refiected in 
the excitation spectrum. 

A. Effective Hamiltonian 

In the center-of-mass frame of the system BEC occurs in the fc = state. We therefore decompose the operators 
that appear in Eq. ( p^ ) into the fe = components and the k ^ ones. The first term on the right-hand side (rhs) 
of Eq. ( 2^ is rewritten as 

■ piP^ PoPo pip'' : + ^ : plPk : • (132) 

k k^O k^O 

If we ignore the terms that do not include the fe = components, we may approximate pk^o as 

Pk^o ^ Y{ai^„Ak.,rn + alfc,™ao,m) =Dk + dI^, (133) 



where 



We thus obtain 



^ : plpk : + ^ : {2DlDk + DkD-u + dIdI^) : ■ (135) 

k k=iO 



Similarly, we may approximate the second term on the rhs of Eq. ( p4[ ) as 

^ : Ffe • Ffc : ~ : F : +2 ^ ^ F • fmn"-k,mO-k,n + X! X! ' fran 

k fe^O rrm k=jtO ijmn 

X (2al^ao^nal^^ak,j + aj,,i4,mafcja-fc," + aojao,naI,ia-fc,m)' (136) 
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where 



1^0, n- 



(137) 



The third term on the rhs of Eq. (g4|) is decomposed into 



where 



(138) 



(139) 



and H.c. denotes Hermitian conjugates of the preceding terms. Substituting Eqs. ( 133 ), (136) and (13S) into Eq. (gj), 
we obtain 



2V 



hV 



+ ^ E E ■ fmnaLrrAk,n + |^ E E ' f ^nni'^4),t^0,nal.^ak,j + aL4,mafc J ^-fe," + «0 jao,naL«-fc,m) 



ky^O inn 



fc^O ijmn 



+^ I ^+ E E(-i)'"'^'^'™«-'^-" + H-c- 

fc^O m 



2C2 



(140) 



k^O mn 



In the Bogohubov approximation we replace operators ao,m by c-numbers VN^^Cm, where Cm's denote the complex 
mean-field amplitudes introduced in Sec. [II and N^^'~^ is the number of condensate bosons. Since N^^'~^ is smaller 
than the total number of bosons N due to the interparticle interactions, we take into account the conservation of the 
total number of bosons through the relation 



TV 



BEC 



E E ^ ^■ 

k^O m 



(141) 



Then Eq. (k3) becomes 



= -pN{f.) + E E [^-^ - ^'("^ - 



(142) 



k^O m 



and Eq. (140) becomes 
V ~N 



'icon , 



K-^-)nEE^'^-^ 



k^O m 



^1'' E E<-^> ■ /mn"l,m«fc,n + ^ E E ' /mn (^C CnaLm"fc J + CCa/cjA-fc^n + OCna^«-/=,m) 



fc^O mn 



/c:^0 ijmn 



C2n 



fc^O m 



H.C. 



2c')n 



EE(-l)"^"CCn4.-m«^,-n. 



(143) 



where n = N/V, and the definitions of (fz), (/), and are given in Eqs. (|l])-(|33|). Equations (|l4^ ) and ( p3|) 

constitute our basic Hamiltonian in the following discussions. We use this Hamiltonian to examine low-lying excitation 
spectra for each phase. 
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B. Excitation spectrum of a ferromagnetic BEC 



Let us first examine the excit ation spec trum of a ferromagnetic phase in which BEC occurs in the m = 2 state. 
Then Cm = 2Sjn,2 and Eqs. (|142D and (|143D become 



Ho = -2pN + J2 [fife + (2 - m)p] nu,m, 



and 



(144) 



, (145) 



respectively, where 54 = co + 4ci . The total Hamiltonian is therefore given by 

1 



-QiuN - 2pN 



-E 

fc^O 



1 



(cfc + 54"-) "-fc,2 + 2.94"(afe,2a-fc,2 + al,2«-fc,2) 



+ E [('^'^ + -P)"fc,i + (efe + 2p - 4cin) fik^Q 

fe^O 

+ (ffc + 3p - 6cin) + (efc + ^p- 8cin + 2c2n/5) nfc,-2] 



2,fc 



254 nck, 



The second line gives the Bogoliubov spectrum 
while other terms give single-particle spectra: 

Em 



For the Bogoliubov excitation energy to be positive, we must have 



(146) 



(147) 



l,k = Cfc - 






(148) 




^2p- 


4cin 


(149) 


= ^k- 


l-3p- 


6cin 


(150) 


2M — f^k - 


|-4p- 


Scirt + 2c2n/5 


(151) 



54 



-04 > 0. 



(152) 



That is, the s-wave scattering length for the total spin-4 channel must be positive. This conditio n is t he same as that 
required for the ferromagnetic mean field to be stable, that is, the first term on the rhs of Eq. ( |146[ ) being positive. 
For the single-particle excitation energies to be positive, we must have p > 2cin and p > (2ci — C2/10)n. These 
conditions are the same as those in (^l|) for which the ferromagnetic phase is the lowest-energy mean field (note that 
p ~ 2p/n). 

We note that the Bogoliubov spectrum (147) is independent of applied magnetic field and remains massless in its 
presence. This Goldstone mode is a conseq uence of the global U(l) gauge invariance due to the conservation of the 
total number of bosons, as discussed in Sec. [II A. 



C. Excitation spectrum of an antiferromagnetic BEC 

Let us next examine the excitation spectrum of an antiferromagnetic phase in which the order parameter is given 
by (^). Making the replacements 



C2 = e*Wi + ^, C-2 = e»^-y^-^, Co = C±i=0, (153) 
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and substituting Eq. (153) into Eq. (143), together with Eq. (142), we obtain the total Hamiltonian of an antiferro- 
magnetic BEC: 



i?^^ = ^(co + C2/5)n7V+ i(ci - C2/2Q)nN{f,f-pN{l) 



E 

fc^O 



2c2n 



ek+piif.) - 2) + Con|C2p + Cin(4|C2p + 2(A) - (L)^) + ^(IC-2p - 2|C2C-2n 



fife, 2 



efe +p((A) + 2) + con|C-2p + cin(4|C-2p - 2(/.) - {^f) + ^(|C2p - 2|C2C-2n 



E 

fe^O 



j54'^(C2^afe,2a-fe,2 + C* 2afe,-2a-fe,-2) + (cq - 4ci + 2c2/5)n(C2Cl2«fe,2a-fe -2 + C2C* 2afc,2^fc -2) + H 
efe +P((A) - 1) + cin(2|C2p + (/.) - (A)^) - ^|C2C-2p nfe,i 



efe +p((/.) + 1) + ci7i(2|C-2p - (A) - (/.)') - ^|C2C-2p nfe.^i 

5 



+2(ci - C2/5)n(C2C*2afe.i«-fe -1 + H.c.)} 

2 4^-- 



E 

fc^O 



ek+p(/.)-Cin(/,)^-^|C2C-2p 



"-fc,0 + — ^(C2C-2«fe,0«-fe.O + H.C. 

5 



(154) 



This Hamiltonian may be simplified using the relation p ~ {fz){ci — C2/20)ri, giving 



^^^ = ^(co + C2/5)nJV-ip7V(/,) 



+ ^ {(Cfc + gin\C,2\^)flk,2 + {(k + 5'4n|C-2|^)fifc,-2 



+ [^(C2^afe,2a-fe,2 + C* l^fc -2a-fe -2) + (cq - 4ci + 2c2/5)7i(C2*C* 2afe,2a-fe -2 + C2C-2aL2"fc,-2) + H.c. | 



E 

fc/O 



efe + (ci - C2/5)n + i^{ci + C2/10)n(/^) 



fife,i 



Cfc + (ci - C2/5)n - -(ci + C2/10)n(/2) 



fife,-i 



+2(ci - C2/5)n(C2C*2afc4«-fe -1 + H.c.)} 
^ |(efc - C2n/5)nkfl + ^((2 C* 2afc,oa-fc.o + H.c.)| 



(155) 



fc^O 



This result shows that the eigenmodes are classified into three categories: the m = Q mode, the coupled m = ±1 
modes, and the coupled m = ±2 modes. Below we analyze each of them. 



1. The spin-0 quasiparticles 



The Hamiltonian (155) shows that the m = mode is decoupled from other modes even in the presence of magnetic 
field. This part of the Hamiltonian can readily be diagonalized to give 



Ek,^Jel+'-\f^ek+{^y{h)^. (156) 



5 " V 10 

We note that the spectrum (|156[) becomes massive in the presence of magnetic field. 



2. The spin-1 quasiparticles 



The rn = 1 and m = —1 modes are coupled in the Hamiltonian (155), and the eigenenergies can be obtained by 
diagonalizing the following Hamiltonian: 
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^ = X! i'^ksfik,! + Cfc -in-fe -1 + (5afc,ia_fe _i + H.c.)} , 



(157) 



where S = 2(ci — C2/5)nQi^'L2 and 



efc,±i = efc + (ci - C2/5)n ± -(ci + C2/10)n(/z) 



(158) 



The dispersion relation can be found by writing down the equations of motion for Ofc i and a^j^. and seeking for 
the sohition of the form exjp{^iEk.±it / h) , with the result 



Ek,±^ - ±1 (ci + ^) + ^4 + 2 (ci - I) n6k + ^ (ci - y)' 



(159) 



The excitation spectra ( 159 ) become massive in the presence of magnetic field. The positivity of this energy is 
guaranteed by the conditions ci — C2/2O > p/2n > and C2 < which are required for the antiferromagnetic phase 
to be the lowest-lying mean field (see condition (|72|)). 



The spin-2 quasiparticles 



The m — 2 and m — —2 modes are coupled in the Hamiltonian (155), and the relevant part of the Hamiltonian 
reads 



h-2 = 



efc,2»^fc,2 + efe,-2"fe,-2 + ~^(C2 ak,2a-k,2 + C-2"fe -aO-fe -2 + H.c.) 



where 



-|-(co - 4ci + 2c2/5)n(C2C* 2afc,2a-fc -2 + C2C*2«L2"fc -2 + H.c.)| 



efc,2 = Efe + .94"-|C2p, efc-2 = efc +g4»T-|C-2| 



(160) 



(161) 



By unitary transformations a±k,2 — > ci±k,2£ and a±fe.-2 ^ fi±fc,-2e*"^~^ , Eq. (160) reduces to 



E 



^2^ = y efe,2"fc,2 + ek-2nk-2 + ^{aak,2^-k,2 + /3afc -2a-fe -2 + H.c.) + j{ak,2a-k-2 + aL,2«fe -2 + H.c.) , (162) 



where a = g4n|C2p: /? = 94n\C-2\'^, and 7 = (cg — 4ci + 2c2/5)n|C2C-2|- This Hamiltonian can be diagonalized by 
writing down the Heisenberg equations of motion for 0^,2 + c^^fc 2 ^^'^ o,k.-2 + a^fc _2 ^ 



{iK)^^{ak,2 + alfc.2) = <^k{£k + 2a)(dk,2 + a^f,.^) + 2^ek{ak-2 + a\^_2), 
{ihY-^{ak-2 + a-k -2) = + /5)(afe,-2 + a^-k-2) + 27efe(afe,2 + a^fc,2)- 

By assuming that ak,±2 + Q-L^ _|_2 oc exp(— ii?^^2^/^)i obtain the following dispersion relations: 

(-^fc,±2)^ = Efc Efc + 94" ± 94" 



At 



, ^ ( C2 

^-^1^^-20 



n 2 



(163) 



(164) 



The positivity of this energy is met if the conditions ci — C2/2O > and cq + C2/5 > are satisfied. The former 
condition is met whenever the antiferromagnetic phase is the lowest-energy state (see (|7^)), while the latter condition 
is required for the antiferromagnetic phase to be mechanically stable, that is, the first term on the rhs of Eq. (155) 
is positive. We note that the dispersion relations (164) are massless even in the presence of the magnetic field. They 
are the Goldstone modes associated with the U(l) gauge symmetry and the relative gauge symmetry (the rotational 
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sym metry about the direction of the apphed magnetic field) that are manifest in the mean-field solution discussed in 



Sec. [II B 



At zero magnetic field, Eqs. (156), ( |15S| ) and (164) reduces to 

«o)' - efe [efc - (2c2/5)n] , 



iEk!±i)^ = [cfc + 2(ci - C2/5)n 



AF ^2 _ j (k [cfc + 2(co + C2/5)n] , 



ek [efc + 8(ci -C2/20)7i] 



(165) 
(166) 

(167) 



implying that all the five excitations are Goldstone modes. This reflects the fact that in the absence of the magnetic 
field the ground state is degenerate with respect to five continuous variables (see Eq. (50)). 



D. Excitation spectrum of a cyclic BEC 

We consider the case of Eq. (^J), namely. 



C±2 = ^(l±a)e'^±% C±i=0, Co = y^^e'^°, (168) 



where a = {fz)/'2 and 02 + 0-2 — 20o = Because (s_) = and (/±) — 0, the interaction Hamiltonian (143) 
reduces to 

k^O m k^O k^O m 

+ 1^ Y Y fijfmniai^tao^nal^ak.j + al^al ,„ak^ja^k,n + H.c.) 

fc^O ijmn 
fe^O mn 



Substituting Eq. (168) into this, performing unitary tra nsfor mations a±k.m ~^ a±fc,me"^™, where Cm = ICmle**^™ and 
02 + 0-2 — 200 = TT, and combining the result with Eq. (142), we obtain 



+ ^ i ^(a + /?)(|C2pafe,2a-fc,2 + |C-2pafc -2a-fc,-2) + ^|Copafc,oa-fc,o + (a - /3)|C2C-2|afc,2a-fc _2 

-|-(a -13 + 7)IC2C-2|aI_2^fe -2 + a|Co|(|C2|afe,2a-fc,o + |C-2|afe,oa-fc -2) 
+ ICo|(a|C2| - 7lC-2|)aL2afe,o + ICo|(a|C-2| - 7lC2|)4,o«fc,-2 + H.c.} 

+ ^ X! {2Co^afc,ia-fc -1 + V6Co(C2afe,ia-fc,i + C2afc -lO-fe-i) + ^{CKo + CoC-2)aI ^lOfca + H.c.j , (170) 

k^O 

where a = cqu, (3 = 4cin, 7 = 2c2n/5, and 

Afc,±2 = efe + (a + /3)|C±2p + 7lCT2|', 
^fc,±i = efe + f(2|C±2p + 3|Con, 
Ak,o^eu + {a + ^)\Co\\ (171) 



It can be seen from the Hamiltonian (170) that there are two separate sets of coupled modes, that is, the m = ±1 
modes and the m = 0, ±2 modes. 
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1. The m = ±1 coupled modes 
The equations of motion governing the m — ±1 coupled modes are given by 



i'h—CLk-i = Ak-iSik-i + L'_2a^fc _i + -BaLfc^i + C*ak.i 



ih—a 



Ak-ia_k -1 — D*_2ak,-i - B*ak,i ~ Ca'_f^^ 



where B = I3Q/2, C = V6/3(C2Co + CoC-2)/4, D±2 = V6/3CoC±2/2. The eigenenergies of these modes are given by 



(172) 



8 



/3 



ffc H 



16 



64 



(173) 



These excitation ene rgies are always positive semidefinite and massive in the presence of magnetic field. 

We note that Eq. (175) has one gapless (but not massless) mode in the presence of external magnetic field. Taking 
the limit |fe| ^ of Eq. (173), we obtain 



(/.) 2 



Sk + ^ifz) and 



(/.)- 



Efc. 



(174) 



However, both of these modes become massless in the absence of magnetic field. In fact, Eq. ( |173| ) then reduces to 

E^^ = Vefc(efc + 4cin). (175) 



2. The m — ±2, coupled modes 



The equations of motion governing the m = ±2, coupled modes are given by 



ih-X+ = {A2 ~{a + l3)Cl)X- - 7lC-2Co|r_ + 7Co^-/2, 

= {A2 + [a + P)Cl)X+ + (2aC2 - lC^2)CoY+ + {a - /3 + j/2)C^Z+, 
ih^Y+ = -7|C-2Co|^- + (Ao - aCl)Y^ - 7lC2Co|^-, 



dt 
d 



(2aC2 - 7C-2)Co^+ + (^0 + aCo)>+ + (2aC2 - 7C2)Co^+, 



in-Z+ = 7Co^-/2 - 7lC-2Co|r- + (A_2 - (a + f3)C-2)Z^. 
ihj^Z, = (a - /3 + 7/2)Co^+ + (2aC2 - 7C2)Co^V + (^-2 + (a + P))Z+, 



where X-^. = 2 i 2; ^± — '^fc.o i a^^ g, and Z^ — 0,^.-2 i a^^. _2- Substituting the expressions for X_, y_ and 
Z_ into those for X+,y+ and Z+, we obtain the equations of motion for the latter set, which reduces to the cubic 
equation and therefore can be solved analytically. The result is given by 



±2.0 ^ 



4+(^)2 / 
Sk + a H z P ± W " 



l--(/,)2)a/3 + 



(4 + (/. 



64 



-/32 



(176) 
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where we recall that a = Cqti, (3 = Acin, 7 = 2c2n/5. The second solutions in Eq. (176) are always positive semidefinite. 
The positivity of the first solution is guaranteed by the condition C2 > 0. We note that the first solution is massive and 
independent of the applied magnetic field and that the second solutions remain massless in the presence of external 
magnetic field. The latter is a consequence of the fact t hat the mean-field solution is degenerate with respect to at 

least two continuous variables, as discussed in Sec. IllC.^ 

In the absence of external magnetic field, the results ( |f 76| ) reduce to 



'^^±2,0 



Cfc + 7: \/efc(efc + 2a), \/ ek{ek + /?)• 



(177) 



It can be shown from the general analytic solutions that these results are valid up to the first order in magnetization 

While we have been unable to complete the analysis of the cyclic phase except for the case of Eq. (|l6§| ) , we would 
like to point out that the excitation spectrum of t he cy clic phase always includes the first solution in Eq. ( |176| ) , even 
when the mean- field solution is not given by Eq. (f6S). This can be seen directly by writing down the equation of 



motion for afe„ 
this equation. 



and 



and the corresponding eigenvalue equation. It can then be seen that -f- 7 is a solution to 



VI. CONCLUSIONS 



In this paper we have studied quantum spin correlations and magnetic response of spin-2 Bose-Einstein condensates 
(BECs) in a mesoscopic regime, and low-lying excitation spectra of each phase of spin-2 BECs in the thermodynamic 
regime. 

The ground states of spin-2 BECs have three distinct phases: ferromagnetic (FM), antiferromagnetic (AF) and 
cyclic (C) phases. The former two phases appear also in spin-I BECs, while the last phase is unique to spin-2 BECs. 
The building block of the AF phase is spin-singlet pairs and that of the C phase is spin-singlet trios. These many-body 
features usually elude mean-field treatments that are based on the order parameter derived from the single-particle 
density matrix. 

Bose symmetry restricts possible building blocks of spin-2 BECs. This can be summarized in terms of the annihi- 
lation operator A^J^^ of n-bosons having total spin /. The fundamental building block is not unique, but one minimal 

set is A2 , Aq , ^2 , Aq , and ' . Bose statistics does not allow units such as Al ' and A\ . The unit ' is 
required to represent a state with odd values of the total spin. 

We have investigated quantum spin correlations and magnetic response in the mesoscopic regime. Under the 
assumption that the system is so tightly confined that the spatial degrees of freedom are frozen, we derived the exact 
many-body ground states which are expressed in terms of the minimal set of creation operators A'^''\ ^o^"*^; ^2'^''^! 
A'Q"\andA'^'"\ These pairwise and trio-wise units help us understand the complicated response of the magnetization 
to the applied magnetic field, which stems from the fact that several values of the magnetization cannot be constructed 
from such units and are hence forbidden. In addition to the quantization of the magnetization to discrete values, 
several new features which elude mean-field treatments are found, such as a sudden jump from the minimum to the 
maximum magnetization, and robustness of the minimum-magnetization state against a small increase in the applied 
magnetic field until it starts to show a linear response. The average Zeeman level populations for the AF-phase ground 
states were calculated, showing that m = 0,±1 populations, which stay zero in MFT, vary sensitively to the applied 
magnetic field. 

We have also examined low-lying excitation spectra using the Bogoliubov approximation. The excitation spectra 
of FM and AF phases are similar to those of the spin-I case In the FM phase, the spectrum consists of one 

massless mode (147) reflecting the global gauge invariance and four single-particle modes (148)- (151) whose energy 
gaps are generated b y th e Zeeman shifts as well as mean-field interactions. In the AF phase, the spectrum includes 
two massless modes ( 164 ) due to the global gauge invariance and the rotational symmetry about the spin quantization 
axis. The remaining three are also Bogoliubov modes, but they all become massive in the presence of magnetic field 
due to th e Zeeman shifts. In the C phase, the spectrum have at least two massless modes (the second term in 
Eq. ( |I76D ) by the same reasons as in the AF phase. The spectrum includes one peculiar single-particle mode (the 
first term in Eq. ( |l76| )) whose energy gap depends solely on the spin-dependent interactions and is insensitive to 
the applied magnetic field. In addition, the spectrum has one gapless mode (the second term in Eq. (174)) whose 
mass dep ends only on magnetization (/z) and vanishes at zero magnetic field. The remaining mode (the first term in 
Eq. (174)) is a Bogoliubov mode which becomes massive in the presence of the magnetic field. 

In the present paper we have studied only static properties of spin-2 BEC. With the very rich phenomena that we 
have found here, we may very well expect that much more remains to be revealed in their dynamics. 
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APPENDIX A: THE ORDER PARAMETER OF SPIN-2 BEC 

The order parameterC of a spin-2 BEC has the same structure as that of the d-wave superconductor which was 
examined by Mermin . We here recapitulate as much of it as is relevant to the our theory. The spin part of the 
order parameter, ^E", of a spin-2 BEC is described as a function of the azimuthal angle and the radial angle 0, and 
may be expanded in terms of the spherical harmonics of rank 2 Y^{9, cj)) as 



(Al) 



m=-2 



where Cm obeys the normalization condition (29). The angle dependence of may be expressed in terms of 
components of a three-dimensional unit vector: — {k^, ky, k^) = (sin cos 0, sin sin 0, cos 0) as follows 



Yr 



±2 



15 
32^ 



{kx i '^ky) , — "Fa/ {kx i iky)kz, Y2 



IGtt 



(2fc2 



^x ^y) 



Substituting these into Eq. ( |AlD , we obtain 



(A2) 



* = 



-n^Mn, 



(A3) 



/ C2 + C- 



M 



1 



Co 



V 



-Ci + C-1 



(A4) 



where 

*(C2-C-2) -Cl+C-1 ^ 

-C2-C-2-y|Co -i(Ci + C-i) • 

-z(Ci+C-i) 2^1(0 j 

The order parameter is thus characterized by a 3 x 3 traceless matrix TrM = with unit normalization 

Tr(M*M) = 1. (A5) 
We may exploit the freedom of the gauge invariance to choose the overall phase so that the real part of TrM^ vanishes. 

ReTrM^ = 0. (A6) 



Let the real and imaginary parts of M be X and Y, respectively. It follows from Eqs. (A5) and (A6) that 



TrX^ = TrY-^ = -. 



(A7) 



Because B is traceless, so can be X and Y traceless. X and Y do not commute, so they cannot be diagonalized 
simultaneously. We follow Mermin to take a representation in which X is diagonal. Then the diagonal elements of X 
become 



27r 



(A8) 



The matrix elements of Y are given by 



Ynn — 



sm 



27r 



n sinx, Y23 = Y32 = 



1 



— - sin (5 cos -i/" cos X, 



1 



Y31 = Y13 = — - sin S sin ip cos x, Y12 = Y21 = - cos 6 cos x 



(A9) 
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Substituting Eqs. (A7) and (A9) into M = X + lY and comparing it with Eq. (A4), we obtain the foUowing charac- 
terization of the order parameter. 

C±2 — ^ (cos 9 ± cos 5 cos X + « cos <j) sin x) 
C±i = ^cosxsinJe^^*'^ 

Co — —/= (sin 6* + z sin X sin cf)) . ( AlO) 
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APPENDIX B: CALCULATION OF ZEEMAN-LEVEL POPULATIONS 



1. Derivation of Eq. (129) 



Let us first write {a^-^am) as 



(0|(5_)^s(5+)A's|0) 



(Bl) 



with 10. 



m/ — "-ml 



What we need is thus a general formula for calculating (■i/;|(5_) (5+) Let us first consider 



the decomposition of \^) into a sum of eigenstates for S , such that 

L»/2J 



(B2) 



where \4>k) is an unnormalized simultaneous eigenstate of {S^,Sz} with eigenvalues S{S—1) and S" = [2{n — 2k) + 5\/A, 
respectively. It follows from Eqs. ( pO| ) and ( ^l| ) that = and Nq = n — 2k. Here n is the number of bosons in \^), 
and [a; J denotes the largest integer that is not larger than x. Let us define ujj and fij such that 



and 



where we have used Eq. (87) and defined the coefhcient (a; b) as 

_ 6!(6 + a + 3/2)! 



(a; 6) 



(a + 3/2)! 



(B3) 



(B4) 



(B5) 



By definition, luq = 1 and ^ /ij = I. Substituting Eq. (B3) into Eq. (p3) yields 



L"/2J 



a;,- 



{n~2j-J)fij+ fik 



{n-2k;k) 



^ '^^ (n - 2/fc; fc - j) ' 



(B6) 



or equivalently, 



L"/2J 

E 



(n - 2j; j) '^'^ (n - 2j;j){n -2k;k~ j) ' 

Using this relation recursively, we can calculate fij as a function of {wfc}(fc = j, j + 1, . . . , \ n/2\ ). On the other hand, 



{n ~ 2k; k) 



(B7) 



multiplying (5+) on both sides of Eq. (B2) and taking norms, we obtain 
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(V-IV) 



{n -2k;l + k) 



/£ = 



{n-2k;k) 



(B8) 



Using this formula, we can evaluate Eq. (Bl). If we apply \Tp) = \<p) , n — s, and / — Ng to the relations (|B3| ), (B7), 
and (B8), and noting that = 0, we have o^o = /io = 1, "j^j = fJ-j = 0(j > 0), and 



(</>|(5_)^^(5+)^3|^) 



(s;7Vs). 



(B9) 



For 



we have iS_|0m) = ['^-,a|; 



(</'l0> 

(— l)™a_m|(?!)) and {S-)^\(j)m) = 0. Using these when we apply 



= s + 1, and I = Ns to the relations (B3), (B7), and (B8), we have oji ~ ((/)|al,„a-m|0)/((/)m|?!'m), 



= 0(j > 1), /ii = uji/{s ~ 1; 1), /^o = 1 - Mil and 



= (s + 1;A's)mo H 7 T-TT Mi- (BIO) 



Combining Eqs. (|B1|), (|B9|), and ( |B10| ), wc obtain 



„l + / „^ (aL^a-m)o \ , (s-l;jVs + l) ..f ^ \ mm 

^ (.;iVs) r™""^" " (.-1;1) + (.;7Vs)(.-l;l)2^^-™"-'"^°' ^^^'^ 



where (• • •)o 



|(/))/((/>|(/>). Substituting Eq.(B5), we obtain Eq. (p9|) 



2. Exact forms of Zeeman populations 

Here we give the resuhs of Zeeman populations (aj,ja„) for the states {A^^^'^)'^^ {^1)"^^ (A^^^'^)"^^ (A^^)'^)"33 |vac) with 

1^22 = Oi 1 and 7133 = 0, 1. 

i) ^22 — and 7133 = 0. We t ake \ (j)) = (a2)"^^|vac) and s — ni2. Then, (0202)0 — ?^12 and (^L'^nOo = for m < 2. 
Putting these into the formula (|129|), we obtain 



(4a2> = ni2 + Ns 



3iVs 
2ni2 + 5' 



(B12) 



{a\ai) 



(aL;^a_i) 



2iVs 



ii) ^22 = 1 and 1133 = 0. We take 



2ni2 + 5' 



{al2a-2) = Ns 



3iVs 



2rii2 + 5 

(a|)"i2[%/3(ai)2 - 2V2a2ao]Vac) and s = ma + 2. Then, 
3 



(a|a2)o = ni2 + 1 



(alai)o 



4ni2 + 7' 



(aS^o)o = 1 



4ni2 + 7' 



4ni2 + 7' 



and {al^am)o = for m < 0. Putting these into the formula (129), we obtain 

5A^s / 2iVs .\ 3 



(4^2) = "-12 + Ns 



2ni2 + 9 V2ni2 + 9 / 4ni2 + 7 



(B13) 
(B14) 

(B15) 

(B16) 
(B17) 

(B18) 
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.f. ^ 2iVs , ( 2Ns 



^^^^^ 



2ni2 + 9 V 4ni2 + 7 

iii) 7222 = and 7133 = 1. We take |(/)) — (02)"^^ [(^^i)"^ ^ \/6a2aiao + 2(02)^0-1]^ |vac) and s = ni2 + 3. Then, 

(aLia_i)o = 1 - (B26) 
2ni2 + 5 



and (0^20-2)0 = 0. Putting these into the formula (12£), we obtain 



/-t' \ STVs f 2Ns , 3(7.12 + 3) 

(a,a2) = 7712 + A^s - ^^^^^ - 1^^;^^^ + 1 j (,^^ + 2)(2ni2 + 5) + ^^^'^ 
= ;^r^ + f TTI^^ + 3(-i2 + 4)) , , (B28) 



27712 + 11 V2r7l2 + ll 7 (7712 + 2) (27712 + 5) 

^ ° 27712 + 11 ^2^12 + 11 ; (7112 + 2) (27712 + 5)' 

(alia-i) = . ^^l^^ + f - 3(7712 + 2)) — i- — + 1, (B30) 

27712 + 11 V2"-i2 + ll / ("•12 + 2)(2ni2 + 5) 

f^f , X _ , 6(7712 + 3) ^^ 



,a_2) = A^s ^ 5 + ^ / ' . (B31) 

' ^ 27712 + 11 V (77i2 + 2)(2r7i2 + 5)y' ^ ' 

iv) 7722 = 1 and 7733 = 1. We take \(^) = (4)"'" [V3(ai)5 - 5^202(01)^00 + 4V3(a2)^ai(ao)^ + 2V3(a2)^(ai)2a_i - 
4\/2(a2)=^ aoa-i]^|vac) and s = ni2 + 5. Then, 

(4a2)o ^ ni2 + 3 - , ^^.T'ltL -.^,V (B32) 

(7112 + 3)(277l2 + 7)(477 + 13) 

^ 9(477^2+ 377712 + 83) 

^"1"^^° - (77i2 + 3)(277i2 + 7)(477 + 13)' ^^^^^ 
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ao«o)o = 1 + 7 -itttt; -irrr-. — -rr, (B34) 

[ni2 + 3)(2ni2 + 7)(4n + 13) 



^'^-^-^^^ ^ ' - K.+3)(2n,2 + 7)(4n+13) ' ^^'^^ 
and (0^2^-2)0 — 0. Putting these into the formula ( |129| ), we obtain 

/~t^ N 77Vs f 2Ns , ^\ 30(ni2 + 4)^ 

(a2a2) = n,, + Ns ~ ^^^-^ - (.^^^I^ + V (n,^ + 3)(2.,2 + 7)(4n + 13) + ^^^'^ 



= ^;^I7TT^ + r + - + + + ''^ ) K2+3)(2n,2 + 7)(4.+ 13) ' ^^^^^ 

(4^0) = + f + 1) 18^2 + 3^,2-3) 

^ ° 2ni2 + 15 V2ni2 + 15 7 (ni2 + 3)(2ni2 + 7)(4n + 13) ' ^ ' 

.t . 4iVs , / , 124iVs , 



aLifl-i) - 7^ ^ + 4iVs(ni2 + 8) - - (8nf2 + 49ni2 + 71; 



2ni2 + 15 \^ ' 2ni2 + 15 ' 7 (ni2 + 3)(27ii2 + 7)(4n + 13) 

(B39) 

2ni2 + 15 V ("-12 + 3)(2ni2 + 7)(4n + 13) / 
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